y=[1,3,5];
u1=[1,3,-2];
u2=[5,1,4];

dot(y,u1)
dot(u1,u1)
dot(y,u2)
dot(u2,u2)

y_W=dot(y,u1)/dot(u1,u1)*u1+dot(y,u2)/dot(u2,u2)*u2
y_perp=y-y_W


